Effect of anisotropy on small magnetic clusters 
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The effect of dipolar interaction and local uniaxial anisotropy on the magnetic response of small 
spin clusters where spins are located on the vertices of icosahedron (Ih), cuboctahedron (Oh), 
tetrahedron (Th) and square geometry have been investigated. We consider the ferromagnetic and 
antiferromagnetic spin- 1/2 and spin-1 Heisenberg model with uniaxial anisotropy and dipolar in- 
teraction and apply numerical exact diagonalization technique in order to study the influence of 
frustration and anisotropy on the ground state properties of the spin-clusters. The ground state 
magnetization, spin-spin correlation and several thermodynamic quantities such as entropy and 
specific heat are calculated as a function of temperature and magnetic field. 



I. INTRODUCTION 

Realizing the promising applications in physics, mag- 
netochemistry and biomedicine, molecular magnets have 
recently been the focal point of intense subject of re- 
search. Although these materials appear as macro- 
scopic objects, i.e., crystals or powders, the intramolec- 
ular magnetic interactions are utterly negligible as com- 
pared to the intramolecular interaction. Thus their mag- 
netic properties mainly reflect the ensemble properties 
of small clusters. It appears that in majority of these 
molecules, the localized single particle magnetic moments 
couple antiferromagnetically and the spectrum is rather 
well described by the Heisenberg model with isotropic 
nearest neighbor interaction, sometimes augmented by 
anisotropic terms and dipolar interaction*^. Thus, the 
interest in the Heisenberg model, which is known for a 
long time, is renewed recently by the successful synthesis 
of new magnetic clusters and magnetic molecules. 

Ab initio studies show that often these magnetic sys- 
tems are frustrated due to the competing magnetic inter- 
actions between the individual magnetic moments. The 
effect of finite sizes, quantum fluctuations and frustra- 
tions can have dramatic consequences on the energy spec- 
tra and can even give rise to new phases apart from the 
conventional Neel-like ovdef^H. A great deal of effort 
has been devoted to the theoretical studies on magnetic 
clusters using diff erent the oretical techniques to solve the 
Heisenberg model 21111211211 . 

Using exact diagonalization of the antiferromagnetic 
Heisenberg model, Konstantinidis et a/P have calculated 
the ground state magnetization for a dodecahedron and 
icosahedron symmetry for s = 1/2 and 1 and have ob- 
tained discontinuity in the field-dependent magnetization 
and double peaks in the temperature-dependent specific 
heat arising due to frustrations. Using perturbation the- 
ory, Coffey et al'?- have studied the effect of frustration 
and connectivity on the magnetic properties of a 60-site 
cluster. Schnalle et al^^ have applied an approximation 
of diagonalization scheme to a cuboctahedron for s = 1 
and 3/2 in order to obtain the energy spectra. In ad- 
dition to the magnetic properties, several studies exist 



for the thermodynamic properties of clusters. For exam- 
ple, Honecker et atf^ have calculated several magneto- 
thermal properties such as the magnetic susceptibility, 
specific heat and magnetic cooling rate for a cuboctahe- 
dron with different spin quantum numbers using the anti- 
ferromagnetic Heisenberg model. Besides the exact diag- 
onalization method, several other te chniqu es such as the 
density matrix renormalization group 2 ^* 221 , cluster expan- 
sions^, spin-wave expansions' 112 ^ and quantum Monte 
Carlo techniques^ 211221 can be used to study the magnetic 
systems. However, some of these techniques have draw- 
backs, for example, quantum Monte Carlo technique has 
limitations in describing the systems with geometric frus- 
tration. The advantage of exact diagonalization method 
relative to these approximate methods is that one obtains 
all informations about the whole energy spectra such as 
the degeneracy, the lowest eigenenergies and eigenfunc- 
tions from which the ground state as well as finite tem- 
perature properties can be calculated. 

In the present work we have applied the exact diagonal- 
ization method to calculate the properties of clusters with 
spin- 1/2 and 1. We have studied the magnetic and ther- 
modynamic properties of small clusters with the focus 
on showing the effect of dipolar interaction and uniaxial 
anisotropy on the magnetization behavior in the presence 
of magnetic field, the studies of which are still limited in 
literature^ 211221 . I n addition, the temperature-dependent 
as well as the ground state spin-spin correlation functions 
are calculated for these clusters and are compared with 
respect to the classical case. 

The paper is organized as follows: In section II, we 
describe the theoretical method used for modeling the 
quantum clusters. Section III discusses the results ob- 
tained for 13-atom clusters with spin-1/2, including the 
effect of dipolar interaction. Then, section IV describes 
the findings for 4- atom clusters with spin-1, where the ef- 
fect of uniaxial anisotropies and temperature-dependent 
correlation functions has been discussed, and section V 
discusses the results for spin-1 icosahedron in the pres- 
ence of local uniaxial anisotropies. In section VI the re- 
sults are summarized. 
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II. THEORETICAL METHOD 

From first-principles calculations it turns outP*l that 
the interaction between electrons may be well represented 
by a model Hamiltonian describing a set of interacting 
spins s*j. An important class of such interacting spin 
models consists of spins coupled bilinearly on a finite lat- 
tice. The Hamiltonian of such a system can be expressed 
by Heisenberg Hamiltonian, 

Hp = - y~] Jjj 8j ■ Sj, (1) 

i<j 

where in general the sum runs over all pairs. Si is the 
spin operator on site i having total spin s and z compo- 
nent of the spin si can take values sf — — s, — s + 1, . . . , s; 
Jij is the exchange interaction. This model describes 
the ferromagnetic (antiferromagnetic) Heisenberg model 
when > (Jy < 0). In one dimension and for only 
nearest-neighbor couplings Jy = J, the s—1/2 Heisen- 
berg model has been solved analytically by means of the 
Bethe ansat^^. Unfortunately, the use of the Bethe 
ansatz is quite limited, as this method is only appli- 
cable to models in one dimension. For higher dimen- 
sions, one has to look for approximate methods. When 
the number of spins in the system is small enough, one 
can solve the problem by employing exact diagonaliza- 
tion techniques 38 . A straightforward way to study the 
model Hamiltonian, defined in Eq. numerically is 
simply to obtain the matrix elements of H in a basis 
of \sf, sf, . . . , sf ), with the z-axis taken as quantization 
direction, where n is the total number of spins in the 
system, and then diagonalize the Hamiltonian matrix nu- 
merically. The Hamiltonian matrix can be decomposed 
into block structure with the use of symmetries of the 
model. Since the isotropic Heisenberg model includes 
only the scalar product between the spins, the Hamilto- 
nian is rotationally invariant in spin space, i.e., it com- 
mutes with the square of the total spin of the system, S 2 
and the z component of the total spin, S z . Even though 
it is straightforward to work in a S z subspace, there is 
no simple method to construct symmetry adopted eigen- 
states of S 2 . Construction of symmetry adopted eigen- 
states of S 2 requires more involved calculations^^. Ad- 
ditionally, the Hamiltonian is symmetric under permuta- 
tions of spins that respect the connectivity of our small 
sized cluster, and the model possesses time reversal sym- 
metry in the absence of external magnetic fields. When 
we take into account the symmetries in the system, the S z 
basis states can be projected onto states that transform 
under specific irreducible representation of the symmetry 
group. In this way, the Hamiltonian is block diagonal- 
ized into smaller matrices and the maximum dimension 
required for numerical diagonalization is considerably re- 
duced compared to full Hilbert space size. 

In the presence of an external magnetic field, dipolar 
interaction and anisotropy the Heisenberg Hamiltonian 



([lj is modified to 

H = Wo — B Z S Z + "Hdipolo + "Hani (2) 

where B z is the homogeneous external magnetic field 
defining, without loss of generality, the z-direction. Here 
it may be noted that the factor g/ig is absorbed into B z 
and z component of the the total spin, S z — sf , can 
take values from —S to S in unit steps, where S is the 
maximum total spin of the system. The dipolar term 
^dipoic in Eq. |2| is defined as 

rtdipoic - T^WMb) 2^ UTTs W 

i<j yl 

where fij = is the unit vector along the line 

connecting the two spins or dipoles located on the sites 
i and j, and the sum runs over all pairs. Hani in Eq. ^ 
represents the local uniaxial anisotropy, which is defined 
by 

Hanx = -^A(el-sl) 2 (4) 

i 

where Di are the local uniaxial anisotropy constants and 
el is the unit vector giving the radial direction from the 
central spin. Since the commutators, [Hdipoie, S z ] 7^ 
and [Hiuu,S z ] 7^ 0, S z is no-more a good quantum 
number in presence of dipolar interaction 'Hdipoio or 
anisotropy term "Hani and therefore in presence of dipolar 
or anisotropy term the block diagonalization with respect 
to different S z values is not possible. 

Now we construct the Hamiltonian matrix in terms 
of eigenstates of the total S z operator and express the 
Hamiltonian in terms of the raising and lowering oper- 
ators sf — sf ± is^. When sf and operates on the 
eigenstates of sf , we have 

stM) = v /s(s + l)-sf(sf ±l)|sf±l>. (5) 
For example, for a spin-1/2 particle, 

tf|ti>=0, sf\U) = \ti), 

Si\ti) = \U), 8T\U)=0, 

and for a spin-1 particle, 

s+|t 4 )=0, a+|0,)=v^|t i > ) 
sf\U) = V2\0i), sr|ti) = V^|0 4 >, 
s i -|0 i )=V2|; i ), S r|;,) = , 

where in the latter case {I, 0, f} denote the three possi- 
ble values of sf = —1, 0, 1. In the absence of ^dipoic an d 
"Hani, the z component of the total spin is a conserved 
quantity and we can decompose the Hamiltonian matrix 
into smaller blocks characterizing each values of the to- 
tal spin. For example, in case of a 13-atom cluster with 
s = 1/2, the 2 13 x 2 13 dimensional Hamiltonian matrix 
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FIG. 1. (Color online) Schematic picture of ICO (left) and 
CUBO (right) with labeling of each atomic site. Both struc- 
tures have 12 vertices with one atom at center (not shown). 
The magnetic field is aligned parallel to r*2 for both ICO and 
CUBO. 



with 



is divided into blocks with dimension ( ,) x ( 1 fe 3 ) 
k = 0, ... ,13. 14 such block matrices have to be diag- 
onalized and the largest block matrix has (g) = 1716 
rows. It may be noted that S z being the good quantum 
number, the Zeeman term is not needed to be included 
in numerical diagonalization process and can be included 
later by shifting the eigenvalues by B Z S Z . However, these 
simplification is not possible when dipolar interaction or 
uniaxial anisotropy term is present in the Hamiltonian. 



III. 13-ATOM CLUSTERS WITH s = 1/2 

We have considered two different geometries, namely 
icosahedron (ICO) and cuboctahedron (CUBO) for the 
investigation of 13-atom clusters. The characteristic fea- 
ture of t he ic osahedron, which has a connectivity like 
fullerenesP^, is that it possesses 12 vertices, 20 triangu- 
lar faces and 30 edges. It is categorized in the symmetry 
group of Ih, which is the point symmetry group with 120 
operations^. On the other hand, a cuboctahedron has 
12 vertices with 8 triangular and 6 square faces and 24 
identical edges, and belongs to the symmetry group Oh- 
Schematic pictures of the ICO and CUBO geometry are 
shown in the left and right panel of Fig. [T] respectively. 
Both geometries possess 12 vertices on the surface shell 
and one atom at the center, and can be transformed into 
each other via a Mackay transformational Though the 
number of nearest neighbors for the center atom are same 
(12) for both geometries, each of the surface atoms for 
both cases possesses different number of nearest neigh- 
bors, i.e., the ICO has 5 and the CUBO has 4 nearest 
neighbors in the outer shell. In this section we first study 
the ground state properties and thermodynamic quanti- 
ties like entropy and specific heat in the absence of dipo- 
lar and uniaxial anisotropy term and then switch on the 
dipolar interaction to investigate its influence on ground 
state properties. 

In the absence of "Hdipoio and Hani, the total Hamilto- 
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FIG. 2. All 8192 energy eigenvalues in units of \J\ for AFM 
interactions of the 13-atom ICO (left panel) and CUBO (right 
panel) with spin-1/2. The energy levels are shaded according 
to their degeneracy. There exists a ±S Z degeneracy in the 
AFM case for both symmetries. For the FM case, the energy 
spectra are reversed with respect to the AFM spectra which 
fulfills Efm = — Eafm- 



nian for a 13-atom cluster with nearest-neighbor interac- 
tion can be written as 



rtis 



-J Si ■ Sj 

i,j>0 



J' 5>V *i (6) 



i>0 



where s*o is the spin of the center atom, and the first sum 
runs over all nearest-neighbor pairs (ij) in the surface 
shell. J is the exchange coupling between atoms in the 
surface shell and J' is the exchange coupling between 
central and surface spins. 

The energy spectrum for the two clusters namely the 
ICO and the CUBO clusters can be obtained by diagonal- 
izing the above Hamiltonian. In the presence of magnetic 
field the ground state energy is obtained by considering 
the minimum of the energy eigenvalues from each mag- 
netization sector. We have calculated the energy eigen- 
values for the different exchange couplings defined by (i) 
all spins ferromagnetic ( J = J' = 1), (ii) all spins antifer- 
romagnetic (J = J' = —1), (iii) central spin is reversed 
with respect to the ferromagnetic surface ones ( J = 1 and 
J' = —1), (iv) antiferromagnetic surface spins with fer- 
romagnetic central spin (J = — 1 and J' = 1). However, 
we will mainly discuss the ferromagnetic (J = J' = 1) 
and antiferromagnetic (J = J' = —1) cases. Note that 
all energies are measured in units of \J\ with \J\ is fixed 
to the value 1 in this work. 

In Table [I] we present closed form expressions for the 
ground state energies E , degeneracies K , as well as the 
lowest energy gap AEi for the different S z sectors of 
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FIG. 3. (Color online) Variation of magnetization M z = (S z ) as a function of the magnetic field for four different exchange 
interactions as listed in the panels of ICO (left panel) and CUBO (right panel). J and J' are the exchange couplings among 
the surface spins and center-surface spins, respectively. The external magnetic field is measured in units of |J|. 



ICO and CUBO. The exact polynomials are determined 
with the Mathematica routine " Root Approximant"^ us- 
ing high precision arithmetics with up to 400 digits. Due 
to the two-fold degeneracy for the ±S Z sector, where the 
minimum energy for each positive S z sector has the same 
value as that of the corresponding negative S z sector, we 
have listed the results as function of \S Z \ only. 

The whole eigenvalue spectrum in the absence of ex- 
ternal magnetic field is depicted in Fig. [2] for the AFM 



TABLE I. Lowest energy eigenvalue Eo, degeneracy Ko and 
lowest energy excitation AEi — Ei — Eo for different S z for 
13-atom AFM ICO and CUBO. Energies are in units of | J\. 



K, 



ICO 



AE[ L 



ciCUBO 



A, 



CUBO 



AEi 



CUBO 



13/2 
11/2 
9/2 
7/2 
5/2 
3/2 
1/2 



21/2 
4 

2-VE/2 
-1.83^ 
-3.96 
-5.420 
-6.28; 



i 



2.382 
0.618 
0.102 
0.0045 
0.022 
0.100 



9 
5/2 


-5/2 
-4.63i[ 
-5.86^ 
-6.06^ 



3 
1 

0.293 
0.339 
0.093 
0.093 



a Zero of x 3 - 5x - 3 = 
b Zero of 64a; 6 + 448a; 5 + 656a; 4 - 1184a; 3 - 3412a; 2 - 2036a; - 
c Zero of 2a; 5 + 16a; 4 + 29a; 3 - 23a; 2 - 61a; -8 = 
d Zero of 4a; 10 + 84a; 9 + 700a; 8 + 2842a; 7 + 4992a; 6 - 1726a; 5 - 

21401a; 4 - 31503a; 3 - 14082a; 2 + 4014a; + 3402 = 
: Zero of 

64. ° 



53 = 



; + 6361 = 



e Zero of 

64a; 6 + 960a; 5 + 4784a; 4 + 7168a; 3 - 8148a; 2 - 19868a; - _ . 
f Zero of 65536a; 16 + 2752512a; 15 + 51707904a; 14 + 571146240a; 13 + 
4089167872a; 12 + 19595452416a; 11 + 61510348800a; 10 + 
109531144192a; 9 + 14047096320a; 8 - 488888621568a; 7 - 
1389656886528a; 6 - 2016792866048a; 5 - 1655926247744a; 4 - 
669806791648a; 3 - 39673588208a; 2 + 46200676992a; + 7484904361 = 
s Zero of 

64a; 21 +3136a; 20 + 70272a; 19 -(-952256a; 18 +8684000a; 17 + 55985680x 1 ' = 
259611872a; 15 + 853909520a; 14 + 1844888624a; 13 + 1761797108a; 12 
3621087792a; 11 -18691236512a; 10 -39464764094a; 9 -49351650308x 
34081746286a; 7 - 3226424608a; 6 + 17175800242a; 5 + 16425687591a;' 



interactions of ICO (left panel) and CUBO (right panel) , 
where the minimum energy eigenvalues for different mag- 
netization of the system are found to be different. More 
over the energy gaps between the minimum energy eigen- 
values of the consecutive S z sectors differ for both ICO 
and CUBO. This observation identifies the influence of 
symmetry on the nature of eigenvalue spectrum of the 
system and also explains the nature of the variation of 
the magnetization with respect to the external magnetic 
field for the two clusters. For the FM interaction, we ob- 
tain degenerate minimum energies. This occurs because 
of the fact that the Hamiltonian has spin rotational in- 
variance and as a result, turning the total spin in another 
direction does not change the energy of the system. 

The left panel of Fig. [3] shows the variation of mag- 
netization in the unit of g/is as a function of external 
magnetic field for the four cases of interactions of the 
ICO (mentioned above). In the presence of an external 
magnetic field, the minimum energy configuration for the 
AFM interaction (J = J' = —1) gives rise to plateaus, 
which have been marked by the solid red curve in the 
left panel of Fig. [3j The appearance of different sizes of 
plateaus is related to the inequivalent energy gaps be- 
tween the minimum energy values of consecutive S z sec- 
tors. On the other hand for FM interaction the ground 
state energy lies in the S z = 13/2 sector for all positive 
values of B z and in the S z = —13/2 sector for all nega- 
tive values of B z . Thus for the ferromagnetic interaction 
(J = J'), irrespective of the values of magnetic field, 
\M Z \ = 13/2 (see the black dashed line in Fig. [3}. The 
right panel of Fig. [3] shows the variation of magnetization 
as a function of magnetic field for the four different set 
of exchange couplings in the case of CUBO geometry. As 
observed in the case of ICO, a similar behavior for the 
variation of magnetization with respect to the external 
magnetic field is observed for the FM interaction. How- 
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FIG. 4. (Color online) Variation of ground state magnetization Mf = (sf} as a function of magnetic field (in units of | J|) for 
the AFM case of 13-atom ICO at several values of the reduced dipole coupling strengths uj. The dark blue lines in all plots 
show the field-dependence of magnetization for the center spin («§)■ The orange and light gray lines show the same quantity 
for the top/bottom atoms and remaining 10 atoms on the surface, respectively. The top/bottom spins are strongly affected by 
dipolar interactions. 



ever, for the AFM interactions of CUBO, the plateaus 
appearing in the magnetization have different sizes com- 
pared to the ICO, which can be noted from the solid red 
curve in the right panel of Fig. [3j The differences in re- 
sults is the consequence of the differences in structural 
symmetries of the two clusters. 

Now we shall study the effect of dipolar interaction 
on the magnetization of the 13 atom ICO with s = 1/2 
in presence of magnetic field. Since dipolar interaction 
breaks the isotropy of the system, dipole-dipole interac- 
tion may be an important source of the observed mag- 
netic anisotropy of various magnetic materials^. More 
over, as dipole-dipole coupling depends only on known 
physical constants and inverse cube of the interatomic 
distances, understanding the role of dipolar interaction 
on different properties of the molecule will be useful in 
the studies of molecular structures. For our studies we 
shall consider the Hamiltonian given in Eq. ^. Using 
the exact diagonalization technique we calculate the mag- 
netic properties of the above mentioned system. Figure [4] 



TABLE II. Ground state expectation values of center and 
surface spins at uj = for the AFM case of ICO and the 
corresponding S z . Note that for S z — 13/2 the central spin 
is oriented parallel and Eq. |8| does not hold. 



S z 



(so) 



13/2 
11/2 
9/2 
7/2 
5/2 
3/2 
1/2 



+0.5000 
-0.4231 
-0.4091 
-0.3889 
-0.3571 
-0.3000 
-0.1667 



0.5000 
0.4936 
0.4091 
0.3241 
0.2381 
0.1500 
0.0556 



shows the variation of magnetization as a function of the 
magnetic field for different reduced dipolar interaction 
strengths, 



/^o (ffMB) 2 
47r|r 0l |3|J|' 



(7) 



w = 0, 0.025, 0.05 and 0.1 for the AFM case, where 
\foi\ denotes the shell radius. For uj = we find a re- 
versed central spin sq with negative hysteresis, as long 
as B Z /\J'\ < 13/2. For larger fields the central spin flips 
into field direction. Table |TT] lists the values of center Sq 
and surface spin sf magnetizations for uj = at different 
values of S z , given by 



S* 



s z 



2{S Z + 1)' 



12 



(8) 



At finite values of uj, the magnetization of the sur- 
face atoms (the light gray curve in Fig. [4| behave dif- 
ferently depending on their position which indicates that 
the dipolar interaction has a strong impact on the magne- 
tization of these spins. The spins of the top and bottom 
atoms of the cluster (see left panel of Fig. [I]) are strongly 
modified (orange curves) compared to the other surface 
spins (light gray) even at very small values of uj, while the 
magnetization of the center atom is nearly unaffected by 
the change in ui values. On the other hand, for the FM 
case, uj has no influence on the magnetization of center 
or surface spins. 

In order to understand the impact of frustration on 
the spin configuration of the clusters with different ge- 
ometries, the spin-spin correlation function for ICO and 
CUBO are also calculated. These are directly connected 
to the magnetic structure factor by a Fourier transforma- 
tion, which in principle can be measured experimentally 
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FIG. 5. (Color online) Ground state correlation functions for n (with 1 < n < 4) shell neighbors of ICO and CUBO for the 
AFM case with spin-1/2 and spin-1. The center and the zero index in the abscissa indicates the correlation functions for the 
center atom and from center to the atom on surface shell, respectively. The dashed lines indicate the classical limit for spin-oo 
ICO and CUBO. For the FM case, the correlation functions for ICO and CUBO with spin-1/2 possess same magnitude (0.25) 
for all neighbors. 



by e.g. neutron scattering techniques. Nevertheless, we 
are not aware of such experiments on clusters. The cor- 
relation function at finite temperature can be defined as 



where, (3 — 1/T with the temperature measured in the 
units of Boltzmann constant. However, at T — the 
correlation functions are calculated from the eigenvectors 
obtained from the exact diagonalization of the Hamil- 
tonian in Eq. ([6| directly. The distance dependence 
of the correlation function at zero temperature are 
plotted in Fig. [5] The correlations for the FM case 
are found to be the same for both geometries, whereas 
for AFM interactions, we obtain different correlation 
functions for ICO and CUBO, suggesting the existence 
of frustration in the system. However, the ICO seems 
to be less frustrated with respect to the CUBO, as a 

regular H 1 — oscillation is found for the ICO, while 

the CUBO exhibits an irregular and smaller correlations 
in the 3rd and 4th neighbor shell. The ground state 
correlation functions for the classical limit s— >oo are 
also calculated for these clusters. While comparing 
the classical and quantum spin correlation functions it 
reveals that the ground state correlation functions for 
the ICO with s=l/2 shows a similar qualitative trend 
compared to the classical case while for the CUBO 
the correlation functions show large deviations in the 
third and fourth neighbors relative to the corresponding 
classical case. This trend indicates that the CUBO has 
stronger effect of frustration compared to the ICO. 

In addition, thermodynamic quantities such as entropy 
S and specific heat C are calculated in the absence of 
dipolar or uniaxial terms in Hamiltonian for AFM in- 



teractions in the ICO and CUBO as a function of mag- 
netic field at several temperatures, which are shown in 
Fig. [6j Sharp peaks at low temperature are observed for 
the AFM case as the magnetic field is changed. This 
is due to the fact that the thermal fluctuation is en- 
hanced at those magnetic fields where level crossing oc- 
curs. With increasing temperature, however, a larger 
number of states from each S z sector contributes to the 
thermodynamics, thereby smearing the peaks of the en- 
tropy S is observed. Similar explanation can be given for 
the behavior of the specific heat with respect to the mag- 
netic field at various temperatures. For the FM case (not 
shown), however, only the maximum S z block matrix has 
the lowest energy for all magnetic fields. In other words, 
all eigenvalues are simply scaled with magnetic field and 
thus trivial features observed in the same thermodynamic 
quantities like S and C and therefore they are not plot- 
ted in the present paper. The thermodynamic observ- 
ables for the antiferromagnetic interactions as a function 
of temperature for both cluster geometries are shown in 
Fig. [7| The peak in the specific heat curve as a function 
of temperature at T ss 1 marks the classical excitations in 
the system. However, both systems also have pure quan- 
tum excitations from the low lying energy levels at lower 
temperatures, visible as additional peaks in the specific 
heat and plateaus in the entropy. For example, in the 
case of ICO at B z = 2, these excitations are at a very 
low temperature T « 10~ 3 and stem from the very small 
energy gap AEi = 0.0045 in the S z — 5/2 sector, see 
Tab. [I] A similar behavior is observed at B z = 1, where 
the maximum at T « 10~ 2 comes from the small energy 
gap A£i = 0.022 in the S z = 3/2 sector. 

Figure [8] shows the variation of magnetization as a 
function of the magnetic field at different temperatures 
for AFM (left) and FM (right) interactions, respectively. 
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FIG. 6. (Color online) Variation of thermodynamic entities as a function of external magnetic field for the AFM case of 
13-atom ICO (left panel) and CUBO (right panel) with spin-1/2. The top and bottom panels show the variation of entropy S 
and specific heat C with respect to the magnetic field B z , respectively, at several temperatures T. 




FIG. 7. (Color online) Entropy S (top) and specific heat C (bottom) as a function of temperature (in units of |J|) for the 
AFM case of 13-atom ICO (left panels) and CUBO (right panels) with spin-1/2 for different magnetic fields B z . 
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FIG. 8. (Color online) Magnetization as a function of magnetic field measured in units of |J| at various temperatures for 
the AFM (left) and FM (right) cases of 13-atom ICO. With increase in temperature (arrows), the plateaus start to vanish for 
the AFM case of 13-atom ICO. The insets for both AFM and FM cases shows the variation of magnetization with respect to 
temperature for several values of magnetic fields (arrows). 



It shows that for both interactions, the magnetization is 
smeared out with increasing temperature. The insets in 
Fig. [8] (left and right) show the variation of magnetiza- 
tion with respect to the temperature at different mag- 
netic fields, which shows that quantum effects vanish at 
around T « 0.1 in the AFM case and that the total 
magnetization tends to decrease with increasing temper- 
ature. For CUBO, a similar variation of magnetization 
as a function of magnetic field is observed at several tem- 
peratures. 



IV. 4-ATOM CLUSTER WITH SPIN-1 

Now we present results for 4-atom clusters with s = 1 
and uniaxial anisotropies (see Eq. Ml)). It may be noted 





FIG. 9. (Color online) Schematic picture of a planar square 
(left panel) and a tetrahedron (right panel). The double ar- 
rows indicate the radial anisotropy axes for s = 1. For the 
square case, the magnetic field is perpendicular to the plane 
and for the tetrahedron, it is aligned parallel to r\. 



that such anisotropies only give a constant for s = 1/2. 
For a spin-1 system, the total Hamiltonian in the pres- 
ence of local uniaxial anisotropy axes e*j reads 



-B z S z -J2d* (el 



■si) 2 , (10) 



where Di are the local uniaxial anisotropy constants and 
el are the easy axes compatible with the symmetry of 
the systerrPSl. We have previously studied the structural 
and magnetic properties of small transition metal clusters 
with more emphasis on the magnetic anisotropy using 
the density functional theory (DFTjP^, where the ener- 
gies obtained from DFT calculations were fitted by using 
a classical Heisenberg Hamiltonian. The investigations 
presented here can be viewed as a continuation of the 
previous work in the sense that we perform exact diago- 
nalization of a corresponding quantum spin Hamiltonian 
to study the system. 

There has been several studies related to the mag netic 
and thermodynamic properties for spin-1 clusters^^ 
through Heisenberg model. However, studies including 
the effect of local uniaxial anisotropy on several prop- 
erties of clusters is still limited^. In the present work, 
we have studied the influence of radial anisotropy on the 
magnetic properties and temperature-dependent correla- 
tion functions for the spin-1 tetrahedron and square as 
shown in Fig. [9] A regular tetrahedron (symmetry group 
T,i) consists of four triangular faces, whereby the trian- 
gles meet at each vertex and are equilateral. A square is 
a regular quadrilateral with D4 symmetry. 

In the presence of radial anisotropy, the Hamiltonian 
is modified to a form as represented in Eq. (10), where 
Di = D is the anisotropy constant and e?j are the easy 
axes which differs for each spin. In Fig. [9] the anisotropy 
axes (double arrows) pointing into the radial directions 



9 



4.0 



2.0 



0.0 



-2.0 



-4.0 



(a) 

- 






- 


FM 




• D= 0.0 




♦ D=-5.0 




■ D= 5.0 







-10.0 -5.0 0.0 5.0 

Magnetic field B z 



10.0 




-5.0 0.0 5.0 

Magnetic field B z 



10.0 



n 2.0 




-10.0 -5.0 0.0 5.0 

Magnetic field 6 Z 



10.0 



4.0 



n 2.0 



| 00 

o 

e 

CD 
CC 



-2.0 



-4.0 







"(d) 






Jr AFM 




-* 1 • D= 0.0 




| ♦D=-5.0 




■ D= 5.0 









-10.0 -5.0 0.0 5.0 

Magnetic field B z 



10.0 



FIG. 10. (Color online) The effect of anisotropy on magnetization as a function of external magnetic field (measured in units 
of | J|) for the spin-1 tetrahedron (a,b) and square (c,d). The results for FM and AFM interactions are shown in the left and 
right panels, respectively. The black line (circles) shows the variation of magnetization as a function of the magnetic field. The 
curves with squares (brown) and diamonds (magenta) show the variation of the same quantity for positive and negative D. 



are shown for the square and tetrahedron. For the 4-atom 
spin-1 cluster, the Hamiltonian matrix is of dimension 
3 4 x 3 4 . For D — 0, the Hamiltonian matrix can be de- 
composed into 9 block matrices with S z = —4, —3, . . . , 4. 
However, at finite D the block matrix structure is de- 
stroyed as the uniaxial anisotropy term does not com- 
mute with total S z and the whole Hamiltonian matrix 
has to be diagonalized. 

The presence of anisotropy D results in a different 
qualitative behavior of magnetization as a function of 
the external magnetic field as shown in Fig. [lO] for the 
tetrahedron (a-b) and the square (c-d). In the absence 
of magnetic anisotropy (D = 0), we obtain a single step 
in the magnetization for FM exchange interaction and 
9 plateaus for AFM exchange interaction for both clus- 
ters, the tetrahedron and the square. The presence of 
anisotropy leads to the smearing of magnetization with 



respect to magnetic field for FM and AFM interactions 
(see the squares and diamonds) of square and tetrahe- 
dron, since "H an i mixes the eigenstates of "Ho with dif- 
ferent total spin values. In particular, for tetrahedron 
geometry, we observe differences in the magnetization 
for positive and negative values of -D, whereas for the 
square geometry, the change of sign in the anisotropy 
does not affect the magnetization significantly. This 
shows the dependence of anisotropy on the structural 
symmetry, which has been observed earlier for 13-atom 
clusters through Monte Carlo simulations^. 

In addition, we have calculated the temperature- 
dependent correlation functions, as defined in Eq. , 
for the 4-atom clusters with spin-1 in the presence of 
anisotropy. The variation of the nearest-neighbor cor- 
relation functions with respect to temperature has been 



plotted in Fig. 11 for the spin-1 tetrahedron (a-b) and 
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FIG. 11. (Color online) Variation of the spin-spin correlation function with temperature (expressed in units of | J|) for the FM 
and AFM spin-1 tetrahedron (a), (b) and square (c)-(f). Diamonds, open symbols and filled symbols represent the correlation 
functions {si ■ Sj) for zero, negative and positive anisotropy constants D, respectively. 



square (c-f) with different anisotropy constants for the 
FM (left panel) and AFM (right panel) case. For both in- 
teractions, it has been observed that the anisotropy mod- 
ifies the correlation function significantly at low temper- 
atures. The correlations are positive in the FM cases and 
negative in the AFM case. For D > the correlations de- 
crease with temperature as expected. However, for neg- 
ative D the correlations are reduced at low temperatures 
due to the quantum effects, as can be seen most easily 
for the square: Here the classical ground state would be 
a FM/AFM state in the direction perpendicular to the 
plane, as in this case all couplings are satisfied and the 
spin directions are all perpendicular to the anisotropy 
axes as required. However, "H ani introduces spin-flip pro- 
cess in the system and as a consequence ground state 
becomes the linear combination of states with different 
magnetization. Hence a strong reduction in correlations 
is observed in presence of anisotropy term at low tem- 
peratures. This reduction of correlations becomes neg- 
ligible at higher temperatures, and correlations decrease 
as usual. While the nearest neighbor correlations in the 
AFM square (Fig. Ill) are negative, the second nearest 
neighbor correlations are positive, as the square is not 
frustrated. 



TABLE III. Size of block matrix, lowest energy eigenvalues 
-Bo, A2?i, and degeneracies Kq, K%, for different S z for 13- 
atom AFM ICO with spin-1. 



s z \ 


Matrix size 


Eo/\J\ 


K 


AEi/|J| 


K, 


13 


1 


+42 


1 






12 


13 


+29 


1 


4.763932 


3 


11 


91 


+17 


1 


4.763932 


3 


10 


442 


+10.763932 


3 


1.236068 


5 


9 


1651 


+5.034063 


5 


0.147456 


1 


8 


5005 


-0.131753 


4 


0.061759 


3 


7 


12727 


-4.663902 


4 


0.032182 


5 


6 


27742 


-8.608201 


3 


0.035204 


4 


5 


52624 


-11.932667 


4 


0.040381 


4 


4 


87802 


-14.679508 


5 


0.011089 


3 


3 


129844 


-16.920343 


5 


0.007569 


4 


2 


171106 


-18.566489 


3 


0.013245 


4 


1 


201643 


-19.506298 


5 


0.007972 


5 





212941 


-19.839976 


3 


0.034983 


3 



SPIN-1 ICOSAHEDRON 



Finally we present the results for the spin-1 icosahe- 
dron with Hamiltonian described in Eq. ([6]). Now the 
Hamiltonian matrix has 3 13 = 1594323 columns and 
rows. For vanishing anisotropy, it can be decomposed 
into block matrices whose sizes are the t rino mial numbers 
( 13 ) 2 for -13 < k < 13, given in Tab. |TTl| together with 



11 



10.0 



c 
o 

'-*—> 

03 
CD 

c 

O) 
03 




-10.0 



-15.0 -10.0 -5.0 0.0 5.0 10.0 15.0 
Magnetic field B z 



FIG. 12. (Color online) Variation of magnetization M z as a 
function of the external magnetic field B z (measured in units 
of |J|) for different values of the radial uniaxial anisotropy D 
for the spin-1 AFM icosahedron. 



the lowest eigenvalues of each S z block. From these re- 
sults we compute the hysteresis curve shown in Fig. 12 as 
black line. The curve is similar to the spin-1/2 case, since 
for fields below \B Z \ < 6 the magnetization varies with 
nearly constant step size from -11 to 11. The plateaus 
at S z — ±11 mark the saturated outer shell, where only 
the central spin still points antiparallel to the external 
field. At larger fields also Sq aligns with the field in two 
steps, one from (sg) = —1—^0 at B z = 12 and one from 
(sg) = 0— >-l at B z — 13. This behavior already shows 
characteristics of the classical limit, in which the steps 
vanish and the magnetization varies continuously with 
field. Nevertheless, also in this limit we find a plateau 
at M z — 11 with reversed central spin, which rotates 
towards the field direction in the range 11 < B z < 13. 

The ground state spin correlation functions for the 
AFM case of ICO with spin-1 are plotted in Fig. [5] (see 
the filled squares). This shows again similar qualitative 
behavior to the spin-1/2 case and show antiferromag- 
netic order in the AFM case, interestingly now the anti- 
correlation between the third neighbor atoms is larger 
than the nearest neighbor value. This shows that the 
correlation approaches towards the classical limit —1, as 
shown in Fig. [5] (the dashed line). On the other hand, 
for the FM case, the correlation functions for all shell 
neighbors become one for the ICO. 

Additionally, we have calculated the hysteresis for the 
icosahedron with radial uniaxial anisotropy according 
to Eq. (10). For these calculations it was again nec- 
essary to work with the whole Hamiltonian matrix, as 
the anisotropy term does not commute with the inter- 



action term and thus destroys the block structure of %. 
Nevertheless, we could calculate the smallest eigenvalues 
and eigenvectors using a Lanczos scheme. The resulting 
curves are shown in Fig. 12 as magenta (D = —5) and 
brown (D = +5) lines, they are similar to the results for 
the tetrahedron with anisotropy shown in Fig. |10|a,b. 



VI. SUMMARY AND OUTLOOK 

We have employed the exact diagonalization technique 
for small spin clusters with 4 and 13 vertices using the 
quantum-mechanical nearest-neighbor Heisenberg model 
and calculated the full energy spectrum numerically as 
well as the ground state energies of the 13-atom systems 
analytically. The magnetic and thermodynamic proper- 
ties as well as the spin-spin correlation functions are de- 
rived from these results. The ground state magnetization 
shows discontinuities accompanied by a magnetization 
plateau as a function of magnetic field for the antifer- 
romagnetic exchange interaction. These magnetization 
plateaus vanishes for temperatures around T > 0.1|J|. 
The ground state correlations suggest that the icosahe- 
dron is less frustrated than the cuboctahedron, since a 

regular H 1 — oscillation is found for the ICO, while 

the CUBO has irregular correlations with smaller values 
in the 3rd and 4th neighbor shell, see Fig. [5j 

We have shown that the dipolar interaction plays a sig- 
nificant role for the magnetization in case of AFM inter- 
action of the ICO. Our investigations show that dipolar 
interactions have a strong influence on the magnetiza- 
tion of surface atoms in an external magnetic field, while 
the field-dependent magnetization of the center atom re- 
mains nearly unchanged by the dipolar interactions. 

The field dependence of magnetization and tempera- 
ture dependence of correlation function on tetrahedron 
and square for s = 1 indicates that the influence of radial 
anisotropies on the magnetic properties strongly depends 
on the structural symmetry of the cluster. 

Finally, we have investigated the 13-atom icosahedron 
with s = 1, which involves quite large Hamiltonian ma- 



trices of dimensions 3 



13 



1 13 



that cannot be decom- 



posed into smaller block matrices if the local uniaxial 
anisotropy axes are present. Using Lanczos methods 
we calculated the lowest eigenvalues and corresponding 
eigenvectors of these large matrices, determined correla- 
tion functions and hysteresis curves and compared these 
results to the spin-1/2 cases. 

Regarding any comparison of this exact diagonaliza- 
tion calculation with experiments one has to first of all 
note that this requires tiny negative exchange coupling to 
observe the quantum effects (such as steps) in reasonable 
magnetic field. Nevertheless, there exists few examples 
in nature such as Mni2-acetate and Mn4-dimer molecules 
embedded in organic ligands, which fulfills this condition, 
see for instance RefsP^I. Another interesting point would 
be a collection of quantum clusters, showing long-range 
order at low temperatures. For such calculations, exact 
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diagonalization results can be used as a basis. This is left 
for future studies. 
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